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The influence of random surface inliomogeneities on spectral properties of open microresonators 
00 ' is studied both theoretically and experimentally. To solve the equations governing the dynamics 

' of electromagnetic fields the method of eigen-mode separation is applied previously developed with 

, reference to inhomogeneous systems subject to arbitrary external static potential. We prove theo- 

' retically that it is the gradient mechanism of wave-surface scattering which is the highly responsible 

^ ' for non-dissipative loss in the resonator. The influence of side-boundary inhomogeneities on the 

Q , resonator spectrum is shown to be described in terms of effective renormalization of mode wave 

• numbers jointly with azimuth indices in the characteristic equation. To study experimentally the 

' effect of inhomogeneities on the resonator spectrum, the method of modeling in the millimeter wave 

OO . range is applied. As a model object we use dielectric disc resonator (DDR) fitted with external in- 

^SJ ■ homogeneities randomly arranged at its side boundary. Experimental results show good agreement 

with theoretical predictions as regards the predominance of the gradient scattering mechanism. It 
is shown theoretically and confirmed in the experiment that TM oscillations in the DDR are less 
affected by surface inhomogeneities than TE oscillations with the same azimuth indices. The DDR 
model chosen for our study as well as characteristic equations obtained thereupon enable one to 
calculate both the eigen- frequencies and the Q-factors of resonance spectral lines to fairly good 



^ , accuracy. The results of calculations agree well with obtained experimental data. 



X 



PACS numbers: 05.40.-a, 02.50.Fz, 42.60.Da 



I. INTRODUCTION 



O , Nowadays microresonators (disk-, ring-, and spherical-shaped) evoke considerable interest because new possibilities 
have recently opened up to develop these types of resonators in the optical frequency range and to utilize them as 
oscillation systems for optical lasers [1]. When used in lasers, such oscillation systems offer a number of serious 
advantages, among which are low-threshold currents, a high quality of the radiation spectrum, etc. Microresonators 
^ ', manufactured as a dielectric disk whose diameter is large as compared to the wave length of the radiation are nothing 
else but the open quasi-optical dielectric disk resonators which have long been known in the resonator technology for 
their potential to effectively sustain the electromagnetic (EM) field inside the resonator volume. The retention of the 
. field is provided due to the total internal reflection (TIR) from the resonator side boundaries of waves making up 
• ■ resonance oscillations. As a result, EM oscillations of whispering gallery (WG) type arise. As far as their excitation is 
^ , not accompanied by the additional dissipative loss, superhigh quality factors are achieved in the DDRs. Specifically, 
QQ ' in laser systems, in spite of DDR's microscopic dimensions (the disk diameter is normally about a few /im large), the 
, quality factors can reach the order of 10^ and even more 

■ The high quality factors of the DDRs, that are often prepared from doped silicon, are generally provided not only 
due to extremely small loss attainable in this material but also governed by resonator geometry and the perfection of 
the crystal it is made of. With a certain number of inhomogeneities (local and/or non-local, random and/or regular) 
^ in the resonator material, the ray picture of EM fields in the resonator changes dramatically as against its perfectly 
homogeneous counterpart. The inhomogeneities give rise to local violation of TIR conditions and in this way result 
in additional energy loss which is evident in the quality factor drop. From the above said there arises the problem of 
studying the effect produced by random inhomogeneities in the DDRs on their spectral characteristics. 

In the theoretic analysis of the radiation loss of the DDR two types of inhomogeneities are normally distinguished. 
To the first type the inhomogeneities of volume nature belong, which are related to regular or random spatial variations 
of the permittivity in the bulk of the material the DDR is made of. The other class of inhomogeneities includes the so 
called "surface" imperfections normally related to the deviation of the DDR shape from the ideal cylindrical one. The 
influence of bulk random inhomogeneities on the resonance spectrum was previously studied by the present authors in 
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the particular case of cavity resonators filled with randomly distributed dielectric particles [3| . It was found that the 
physical mechanism through which inhomogeneities affect resonator spectrum is basically the intermode scattering. 
In the case of quasi-optical cavity resonator the peculiar feature of this type of scattering is the selective impact of 
inhomogeneities on different resonance lines. The most affected lines appear to be those which are the least separated 
on the frequency axis, whereas solitary lines are subjected to much lesser influence. Owing to such a selective effect of 
random inhomogeneities, the originally dense spectrum of the quasi-optical cavity resonator is considerably rarefied. 

In Refs. [3j U S [3 1 the investigations were undertaken into the impact of surface inhomogeneities on spectral 
properties of open dielectric resonators of cylindrical and spherical shape. The influence of boundary roughness 
upon the resonance lines was described by means of the simplified quasi-geometric approach where EM oscillations 
scattering due to edge inhomogeneities is taken into account by incorporating into the wave equation the fictitious 
polarization currents (PC) randomly distributed in space _8J. In particular, within the framework of the model 
adopted in Ref. [J| the electromagnetic fields close to the resonator side surface were considered as being excited 
by randomly distributed near-surface current sources whose physical parameters were phenomenologically expressed 
through statistical characteristics of boundary asperities. It was precisely the radiation produced by these sources 
that led to the radiation loss of the resonator and, consequently, to the quality factor reduction. 

Such an essentially phenomenological approach to the description of the effect produced by the roughness of open 
resonator boundaries on its spectral properties cannot be reckoned as satisfactory. The volume current method 
suggested in Ref. \§\, which formed the basis for the PC concept, is, to a large extent, rough and approximate. In 
this method, the radiation loss is most frequently calculated using Green function of the Helmholtz equation. The 
particular form of this function is normally chosen proceeding from the prospective solution of wave equation in the 
far wave zone. Meanwhile, the very notion of the far zone is poorly defined for multiple sources located around the 
periphery of the quasi-optical DDR we deal with in this particular work. Specifically, when deriving characteristic 
equations for open DDR eigen-frequencies, one needs to join EM field components exactly at the boundaries of the 
system under consideration. At this point a significant uncertainty can arise because the actual external fields subject 
to local joining with the internal ones in the presence of surface roughness can deviate considerably from the fields 
approximated into the boundary vicinity from large distances, i. e. from the far wave zone. 

To correctly determine the EM fields near the random-inhomogcneous resonator surface the theories specially 
adapted for the description of classic and/or quantum wave scattering by rough interfaces should be applied (see, 
e.g., Refs. [1, Hi El ) ■ These theories are generally applicable to the cases where wave scattering by random rough 
surfaces is in a sense weak. Typically, this implies fiuctuations of system boundaries to be relatively small in height and 
sufficiently smooth, so that the applicability of Rayleigh hypothesis [l^, [l^ was not violated. However, for confined 
systems like the DDR considered in this work the issues pertinent to wave scattering by surface inhomogeneities 
appear to be much more complicated. Firstly, in practice the inhomogeneities are not always small enough, as well as 
sufficiently smooth, therefore the conditions for scattering weakness are easily violated. Moreover, as the resonance ray 
trajectories in high-Q resonators are periodic, the effect of oscillation scattering caused by the boundary roughness 
is "path" -accumulated. Scattering can become strong even though the asperities are small-in- height and smooth. 
This makes the applicability of the above-mentioned theories of rough-surface scattering in the case of high-quality 
resonators highl y q uestionable. 

In Refs. |14l . Il5j . the novel transport theory was developed for waveguide-shaped systems with random rough 
boundaries. In the framework of this theory it was revealed that the wave scattering resulting from fluctuations of 
inter-media surfaces can be efficiently described in terms of two physical mechanisms, namely, the amplitude and the 
gradient ones. For the first mechanism it is just the mean-square height of the asperities that serves as a main guiding 
parameter, whereas for the second one the mean slope of the asperities or, in other words, their sharpness, plays the 
decisive role. Both of these mechanisms contribute additively to the scattering amplitude, but partial probabilities 
pertaining to them may differ essentially. It was shown in [3, [l^ that, in most cases, the role of the gradient 
scattering appears to be prevalent. Yet there have been no experimental confirmations of this fact in the literature 
so far. 

Apart from peculiar problems associated with surface nature of scattering in edge-disordered resonators, some more 
questions arise which need to be resolved when studying the spectra of DDRs subject to surface inhomogeneities. 
Among these questions, particularly, are those related to the vector nature of fields which are excited in these 
essentially non-integrable systems. It is known that even for ideal cylindrical DDRs the electrical- and magnetic- 
type oscillations cannot be separated in the strict sense. They always remain to the certain extent intermixed 16]. 
This fact renders theoretic analysis of such systems spectra rather sophisticated. Spectral properties of DDRs without 
random inhomogeneities were examined in a number of theoretical papers (see, e.g., Refs. Yet, in deriving 

characteristic equations some inexact a priori assumptions where used, which still require both theoretical grounding 
and experimental corroboration. 

In our present study, one of the goals is to investigate theoretically the physical mechanisms responsible for widening 
resonance lines of dielectric microresonators with random surface inhomogeneities. In view of such systems being non- 
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integrable, it is necessary to elaborate an appropriate theoretical model allowing for sufficiently accurate determination 
of both the frequency spectrum and the quality factors of resonance lines. Yet another goal of this study is to 
corroborate experimentally which of the physical mechanisms does play a dominant role in the scattering of EM 
oscillations excited in microresonators with random rough side boundaries. 

From the theoretical viewpoint, the spectrum analysis in our study is carried out in terms of scalar potentials, 
specifically, electrical and magnetic Hertz functions [i6|. We formulate the conditions wherein these potentials make 
independent contributions to EM fields in the resonator, which is equivalent to decoupling the oscillations of TE and 
TM polarization. The Helmholtz equation for Hertz potentials in an irregular-shaped open resonator is equivalent 
to Schrodinger equation for electrons moving in the piecewise continuous space subject to random potential. Owing 
to this one can extend the results obtained in the present study to quantum systems as well, in particular, to open 
quantum dots having random rough boundaries. 

To solve wave equations in the rough-bounded resonator, the eigen-mode separation method is used, previously de- 
veloped with reference to waveguide-like systems subject to arbitrary static potential [lOj HH US] ) ■ Using this method, 
the originally posed statistical problem of determining the fields in three-dimensional DDR with complex, randomly 
rough, side boundary is rigorously reduced to the set of one-dimensional dynamic equations that contain some ef- 
fective random potentials. We show that under the conditions where gradient scattering mechanism is dominant, 
the boundary roughness effect on the resonator spectrum can be described through certain renormalization of mode 
wave numbers and azimuth indices of the Bessel functions in the characteristic equation. The values of renormalized 
wave numbers and mode indices decrease as the asperities get sharpened, which is consistent with a decrease in the 
resonator Q-factor. We are thus led to the conclusion that the observed reduction of resonant line quality factors 
results not from the extra dissipative loss but rather from EM field intermode Rayleigh scattering induced by random 
surface inhomogeneities. The imperfection of the resonator shape results in the local violation of the TIR conditions. 
This leads to additional radiation loss of EM energy and, hence, to a decrease in the level of localization of EM fields 
inside the resonator. 

Since characteristic dimensions of surface inhomogeneities in actual microresonators are always quite small (of the 
order of nanometers), to verify our theoretical findings experimentally we have decided upon the method of simulation 
with macroscopic devices. As a model system we have employed a millimeter wave quasi-optical resonator made of a 
circular tefion disk. WG oscillations of TE and TM types were excited in the disk using a special waveguide antenna. 
The inhomogeneities of the resonator side boundary were made in the form of teflon bracket-bars randomly attached to 
the outside cylindrical surface. Our experimental results have demonstrated excellent qualitative agreement with the 
developed theory as regards spectral lines widening caused by the resonator side boundary roughness. Furthermore, 
the relatively simple model of the DDR field distribution, which was adopted in our study, enabled us to calculate both 
the Q-factors and the frequencies of the resonance lines with quite satisfactory accuracy. The calculations appeared 
to be in fair conformity with our experimental data. 



II. THEORETICAL MODEL AND DERIVATION OF BASIC EQUATIONS 

Consider an open disk resonator as a finite-height cylinder made of the dielectric material with permittivity Eq- 
Plain-parallel end boundaries of the cylinder traverse the central axis (z) at points z± — zLH/2 (see Fig. [T|), the side 
boundary {S) is formed by the generatrix passing parallel to z-axis along closed contour C whose distance from the 
central axis is given by 

r{ip) ^ R + ^{^) , (^g[-7r,7r], (1) 

R is the radius of the averaged (i. e., ideally circular) contour Co containing no random bends. The function S,{ip) will 
be regarded as the Gaussian random process with a zero mean value, (^(v?)) — 0, and binary correlation function 

{am^'))^'j'w{v-v') . (2) 

Here, a is the mean-square height of boundary asperities, W{ip) is the dimensionless function which has a unit maximal 
value at zero argument and falls to parametrically small values at angle distance ipc I (the correlation angle) . The 
angle brackets in Eq. ^ denote statistical averaging over the ensemble of realizations of random function ^{ip). In 
what follows, for the estimation purposes, along with angle parameter ipc we will use another correlation parameter, 
viz., arc correlation length Sc = ^PcR- Both the random function ^(<p) and the regular W{'^) will be thought of as 
periodic with period 27r. 

It is well known that in general the EM field of an open dielectric resonator can be given as a superposition 
of oscillations of electric and magnetic types (TM and TE polarized, respectively) [Ifjj. The vector fields of both 
of these types can be expressed in terms of scalar potentials, U{y) and V{y) (the so-called electric and magnetic 




FIG. 1: (Color online) The sketch of the dielectric cylindrical resonator with a randomly rough 

side wall. Vectors n and no point out the local normal directions to the rough (S) and to the 
smoothed resonator side surface, respectively. 



Debye potentials), which meet the same Helmholtz equation but are subject to different joining conditions at the 
resonator boundaries. In arbitrarily shaped resonators TM and TE oscillations are essentially intermixed, being 
strictly decoupled only in the case of sufficiently symmetric systems such as, e. g., an infinite dielectric cylinder [l^ or 
a dielectric sphere [6*] . Below it will be shown that TM and TE oscillations can also be decoupled (with high accuracy, 
though approximately) in the case of the cylindrical DDR of arbitrary thickness. This fact opens an opportunity to 
extend the conclusions of the present work to electronic microresonance systems, e.g., partially open quantum dots. 

In studying oscillations in the DDR with randomly rough side boundaries we first reduce the problem of wave 
scattering at the boundary to the problem of scattering in the bulk of the resonator of ideal circular form. The 
Helmholtz equation for scalar wave field \E'(r), whose role is played by one of the above-mentioned potentials, after 
rewriting it in cylindrical coordinates and using the conformal coordinate transformation. 



z , 



(3) 



is reduced to the form 



1 d 



1 ^ + K\r,z)-V('^^-V('^^ 



*(r, z) = 



r dr dr dip'^ dz^ 

(the tilde signs over coordinate variables from here on arc omitted). Here, K^{r, z) = k^e{r, z), 

£o +'i'/To , (r, z) e 17 
1 , {r,z)^n ' 



e(r, z) 



(4) 



(5) 



1 /tq is the phenomenological frequency parameter which takes into account dissipative and other uncontrollable losses 
in the system, except for the radiation loss, is the bulk region occupied by the dielectric. Effective potentials V^^'' 
and y^**^ in Eq. ([4]) are the operators whose coordinate representation reads 



VW = 



(92 

oz'^ 



1 



Civ) d 

R(3{ip) dip dip Rf3{ip) 



" 1 d 




r dr 





ld_ d_ 

r dr dr 



(6a) 
(6b) 



Indices "/i" and "s" specifying potentials ^ indicate that the corresponding potential is mainly governed either by 
fluctuations in the asperity height (i.e. by height function ^{(p)) or by fluctuations in the asperity slope (i.e., by 
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slope function ^'{ip)). Such a subdivision of the potentials, with regard to the asymptotic suppression of correlations 
between functions and Civ) over angle intervals that are large as compared to the correlation angle, urges one to 
consider potentials ([5a|) and as corresponding to different physical wave-surface scattering mechanisms ■ 
We will below refer to these mechanisms as the amplitude (or the height, h) and the gradient (or the slope, s) 
scattering mechanisms, respectively. 

Subsequently we will examine the boundary asperities which are sufficiently small in hight so as to meet inequality 

CT < i? . (7) 

In contrast to the widespread belief (see, e.g., Ref. [l]), this docs not necessarily imply that the surface-roughness- 
induced scattering be regarded as weak. We note, for example, that the local value of the potential (|6b|) is estimated 
by the parameter cr/sc, which may take an arbitrary absolute value provided the condition ([7]) is met. The true 
conditions for the scattering to be classified as weak will be provided below, based on the operator technique applied 
to perform the principal calculations. 



III. SEPARATION OF AZIMUTH MODES IN RANDOMLY ROUGH DISK RESONATOR 



The potentials V'^'^^ and V^^'' in Eq. (H]), that account for inhomogeneity of the resonator side-boundary, are defined 
not quite conveniently from the viewpoint of subsequent use of perturbation theories. The inconvenience relates to 
non-zero average value of the potential (j6bp . By separating this average, we can rewrite Eq. ^ as 



r or or OLp'^ Oz^ 



'4j{r,(p,z) = 



(8) 



where two different slope potentials are introduced instead of potential Eq. (j6bp . which have zero mean values, viz 



R 



ld_ 
r dr 
ld_ d_ 

r dr dr 



The parameter S in Eqs. ([8]) and (j9bp . which is defined as 



(9a) 
(9b) 



i?2 



(10) 



serves as a measure for the mean-square slope of boundary asperities. 

Considering that irregularities of the resonator side surface are aligned paraxially, they cannot result in additional 
scattering along axis z except for the partial reflection produced by the initial step change in this direction of the 
permittivity of the medium. This makes it possible to remove from Eq. ([5]) the dependence on the coordinate z with 
the model method similar to the one we use when analyzing the spectrum of homogeneity-free DDR (see Appendix E| . 
Specifically, we assume the dependence of the mode wave function on the coordinate z to have the same form used 
to find s- and a-solutions for the ideal-shaped resonator (for the particular iJ^-symmetric solution this dependence 
is given in Eqs. (jA2[) and (|A3p ). Afterwards, equation ([8]), being divided by the factor of f + S^, is reduced to the 
following form. 



1 a ^ a 



1 



r dr dr (1 + ^^)r^ 



/92 _ _ ~(s2y 

+ Kl{r\z^)-V^''\z^)-V -V 



4'(r,(/j|zg) = 



where 



Kl{r\z^) 



[ee{R-r)+e{r~R)\k'^-kl , 

1 + S2 " 

WKzi) = -Kl{r\z^)[P^{^)-l] 



\z\ <H/2 
\z\ > H/2 ' 



(11) 

(12a) 
(12b) 



The tilde signs over the potentials in Eq. (fTT|) denote their renormalization by the factor of (l -I- ^ . Equation pT|) 
actually describes two-dimensional fields because symbol z^ acts not as the current axial coordinate but stands for 
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some axial index specifying in which domain of the z-axis — either \z\ < H/2 or \z\ > H/2 — the sohition to Eq. ^ 
is sought. 

In the azimuth mode representation, equation (|lip assumes the form 

*«(r|z5) - Un-^{r\z^)^m{r\z^) = . (13) 

ra^n 

Here, 

^nm(f|%) = <j> dip{ip,n\V{r,Lp\z^)\ip,7n) (14) 

~ ~{s2) 

is the matrix element of the entire potential V{r,(p\z^) = '{z^) +V +V , which is taken between eigen- 
fmictions of the azimuth part of the Laplace operator (see Appendix V,i(r|zj) = Unn{T\z^), the index n changes 
to n = n/y/1 + S^. Later we will refer to matrix elements V„ and Unm as the intra- and inter-mode potentials, 
respectively. 

In the general case it is rather difhcult to immediately solve the infinite set of coupled equations p3p . However, the 
solution can be obtained in terms of the operator technique previously developed by the present authors with reference 
to waveguide-type random systems of arbitrary dimensionality 20, _21J and advantageously applied afterwards to the 
analysis of bulk-disordered cavity resonators [3|. The above technique, as applied to the problem touched upon in 
the present paper, is adapted in Appendix [Bl The advantage of the technique is that it can be used to derive precise 
closed equations for wave functions of each of the azimuth modes, namely, 

*„(r|zs)=0. (15) 

Here, along with local intra- mode potential V„, the operator potential 7^ arises, which rigorously allows for the inter- 
mode scattering. The structure of this potential, though well-recognized, is quite complicated to be operated with 
at an arbitrary scattering intensity. Yet, the estimations we provide in the next section substantially simplify the 
T-potential in Eq. (|15p in different limiting cases, and in this way they allow one to obtain the oscillation spectrum 
of a randomly rough DDR in almost all physically sensible situations. 



IV. SPECTRUM OF THE DDR WITH WEAKLY ROUGH SIDE BOUNDARY 



Hereinafter we will refer to the system as weakly rough one if r.m.s. height of its boundary asperities meets 
inequality ([7]). As the potentials in Eg. (|15p are of the operator nature, we will estimate their strength using standard 
definition of the operator norm [23l [2J| . For our purposes the formula 

^2. ((AV^A^)) .... 

(A ) = sup — — — — (16) 

is the most appropriate, in which the parenthesis symbolize the scalar product on the functional space X consisting 
of the class of solutions to Eq. with no random potentials. 

To estimate the potentials in Eq. (jlip we first calculate their azimuth matrix elements. For relatively small-height 
asperities the "height" potential F^'*^ equals approximately V^'^^z^) « —2K'j_{r\z^)£,{(p)/R. Mode matrix elements 
of this potential are 

^(n) is the Fourier transform of the function S,{ip). Matrix elements of the "slope" potentials V and V are 
equal, respectively, 

Ui^^Hr) = - „,,^ in'-m')iin-m)--^ , (18a) 

\/27r(l -I- ^^)R r Or 

1 1 /-I Pi 

^ 2n{l + E^)R^ E (" - - -) [an - m - m) - (c(n - 1)^1 - m))] --r- . (18b) 



7 



A. The efficiency of intra-mode scattering 



We will ignore the role of the potential p7|) for the intra-mode scattering proceeding from the following consid- 
erations. It is evident that under condition ^ the uniform azimuth mode of this potential can result, at most, in 
relatively small ^/R) renormalization of the unperturbed "in-plane energy" K'^{r\z). Yet, in the case we consider 
below even this does not occur. We will regard the asperities as not only being small in height but also as small-scaled 
ones in the sense that the inequality holds 



< R 



< 1 



(19) 

can be put equal 



In this limit the imiform mode of the random process with parametric accuracy (^(0) ~ a^p, 
to zero, since the process is nearly ergodic within the interval of the angle variable change. 

Gradient potential (|18ap is not involved in the intra-mode scattering by its definition, so that only the potential 
(jl8bp should be taken into account, where one must let m — n. We use as trial functions in Eq. (fT6| the solutions 
of Eq. with random potentials equal to zero, i.e., Bessel function Js [i^±(r|zg)r] in the interval < r < i? and 
Hankel function H-^'^ \K ±{r\z^)r\ in the domain i? < r < oo. Such a choice is justified provided the intra-mode 

scattering has a slight impact on the mode energies. Then we arrive at the following estimate for the potential V„ 
norm. 



m\") 



fl + S2)2i?4 



OO 



/— — oo 



(20) 



After averaging Eq. (j20p using correlation equality 

{mC{l))^V2^<^^W{k)Ski , (21) 
which immediately stems from Eq. ([2]) {W is the Fourier transform of correlation function W{ip)), we obtain 



(llVnf) 



(1+S2)2 \R 



R{i + aysi) 



(22) 



It can be readily seen that the factor standing at K'^ in the right-hand side of Eq. (|22p is, in view of Eq. ([7]), small 
as compared to unity. This substantiates the above assumption about the intra-mode scattering weakness. 



B. Comparative estimations of inter- mode potentials 

The inter- mode scattering rate in our resonator is up to the norm of the operator R entering T-matrix (jBlOp . 
1. The "amplitude" inter-mode scattering. 

When estimating the norm of the "height" item R^'*^ in the operator R one is faced with a need to evaluate the 
expression 

(llR^'^f ) - ^E/^'^^ E / '^i^'-i / r2dr,G^r{r,r,)G^P\r,r,){ui':^^{r^)U^^^^^ 

fc g mi,m2#fc Q 

(23) 

(in order not to overload subsequent formulae we omit axial index z^, as this cannot lead to misunderstanding; function 
G'^^''(r, ri) is introduced in Appendix iBjl . The simultaneous presence in Eq. of both the integrals over radial 
coordinate r and the sums over mode indices stems from the definition of the functional space where the operator R 
is effective. 

In view of Eqs. (fT7|) and (|2ip the correlator in the integrand of Eq. (|23p is calculated to 

(Z^l™\(^i)^i™:(^2)) =2y|(^)'x2(,^)^2*(^^)^(^_^^^)^^^^^^ (24) 

By substituting this into Eq. (|23p and assuming the asperity correlation function to have the Gaussian form, viz. 
W{(p) — exp ( — (/?2/2(p2j , we arrive at the following estimate for the height term in the inter-mode scattering operator, 

(llRC^^f ) - (^iCT)2 . (25) 
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The Rayleigh parameter ka, to whose square the right-hand side of Eq. (|25p is proportional, may take small or large 
values, depending upon the relationship between mean-square height of the asperities and the wave length of the 
excited oscillations. According to the value of this parameter we will distinguish between weak {ka <^ 1) and strong 
{ka 1) inter-mode scattering caused by the "height" potential. 
2. The "gradient" inter-mode scattering. 

This type of scattering between different azimuth modes is related to the availability in Eq. (fT3|) of the potentials (jl8|) . 
Their correlators needed to estimate the operator norms of items R'^''^^ and R(''2) in the operator R, are equal, 
respectively. 



Id d \ f 1 d d 



(k - li){h " mi){k - l2){l2 - m2)(e(fc - ~ rn^)C{k - l2)C {h - m^)) . 

(26b) 



X 

il 5^2 — — OO 



Assuming the function ^{f) to be the Gaussian random process, the equality (j26bp can be partially simplified as the 
double sum over li and I2 is reduced to the single one, 

00 oc 

^ . . . = AiraH^.m, E - ^)'(^ - mifW{k - l)W{l - mi) . (27) 

Then the alternate substitution of correlation functions (PB)) into Eq. (|23p in place of the height potential correlator 
results in the following estimation formulas for operators R'^^^ and R^''^^ 

^ll^''ll^)-(rw-4' ^''^^ 

{\\^''''^^)-{jT^)\K^Rf . (28b) 

The collation of norms (PS)) and reveals that the height and the gradient scattering mechanisms, given their 
effect is evaluated as a function of roughness statistical parameters and radial components of the azimuth mode wave 
vectors, can essentially compete against one another. However, the general statement boils down to the fact that 
it is the gradient scattering which is prevalent in the most part of the parameter region. The amplitude scattering 
mechanism can also dominate, but this is the case for quite smooth asperities only and for oscillations of very large 
azimuth indices. One can conclude from Eqs. ((25|) and ((28| that for the amplitude scattering to be dominant it is 
necessary, first, that the average tangent of the asperity angle of slope obey the condition 

1/2 



tane* < — • (29a) 

Secondly, the inequahty must be fulfilled {K±R)^ » {R/sc)'^, which is admissible for oscillations with quite large 
mode indices, i. e., 

R^' 



n > . (29b) 

If only one of the inequalities ((29)) fails to hold (indeed, this is the case for the overwhelming majority of parameters 
from the above-mentioned set), the gradient scattering appears to be dominating. 



C. Whispering gallery modes of the rough-side DDR 



Based on the above estimations, equation that governs the azimuth mode spectrum of the resonance system 
under study can be substantially simplified if one considers the limiting cases of weak and strong scattering. As is 
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seen from Eq. (|22|) . for small r.m.s. height of the asperities (in the sense of inequality ([7])) the intra-mode scattering 
due to potential Vn is weak, thereby making it possible to disregard it in the main approximation. 

As far as the inter-mode scattering is concerned, it is not straightforward to estimate it in the same way as the 
intra-mode one. From Eqs. ([25|) and ([28|) it follows that inter-mode scattering should be classified either as weak or 
strong one depending upon whether the operator R norm is small or large as compared with unity. If the inter-mode 
scattering resulting from both the amplitude and the gradient potentials is thought of as weak one, which occurs 
when two inequalities hold simultaneously 



a Sc 
Sc R 



(30a) 
(30b) 



the potential 7^ in Eq. pS]) can be ignored with parametric accuracy. In this instance the resonator spectrum can 
be obtained from the equation which differs from the initial unperturbed one simply by rcnormalizing both the mode 
index (n n) and the radial wave number {K^^ K±). In the parameter region where both of the inequalities (|30p 
are satisfied this renormalization is extremely small (S <C ^Pc) and can be safely disregarded. 

Conditions ((30)) . that indicate the inter-mode scattering weakness, can be easily violated. By virtue of peculiar 
technology, when manufacturing micro- and nano-sized quantum resonance systems it is normally difficult to obey, 
e.g., inequality (j30bp which corresponds to extreme smoothness of surface inhomogeneities. Moreover, inequality 
pOap is fulfilled within only limited, i. e., long- wavelength, part of the resonator bandwidth. 

If, at least, one of the conditions pO|) is violated, the inter-mode scattering fails to be classified as a weak one in the 
sense that the norm of scattering operator R in the potential (jBlOp becomes large as against the unity. Nevertheless, 
in this case one may as well simplify the T-potential by expanding it in a series in the inverse scattering operator, 
R^^. The expression between the projection operators in Eq. (jBlOp is identically transformed as 



R 



U+U[t-R] ^ = -U -UR-^ (t-R-^^ ^ 



(31) 



Here we use the symbolic notation Q^^^~^ for the full (i. e., including the whole set of azimuth modes) Green operator 
with no intermode potentials. As a result of "coating" Eq. ((3T|) with the projection operators P„ the first term in 
its right-hand side is eliminated, and the operator acting on the wave function in Eq. (|15p assumes the following 
approximate form. 



1 d d j>2 / I N 



2G(^)-i+P„^(^)"iR-ip„«2 



I d d r^T / I N n 
r or or ^ r^ 



(32) 



By comparing the expression in the right-hand side of Eq. (|32p with the one in its left-hand side, with potentials 
Vn(r|zj) and 7^ set to zero, one can notice that the limiting cases of weak and strong (in the sense of the operator R 
norm) inter-mode scattering differ from one another solely by doubling the wave operator in the latter case. Clearly, 
such doubling implies that in strong surface scattering the amplitude of the excited oscillations is twice as small as 
the one in the weak scattering limit. However, it is evident that the change in the common factor multiplying the 
wave operator cannot reveal itself in dispersion relations. 

We thus have demonstrated that the formal algebraic structure of the wave operator remains basically unchanged 
under conditions of weak and strong scattering. In both of these cases the wave operator differs slightly from its unper- 
turbed form. The main difference between dispersion relations for resonators with perfect and rough side boundaries 
is that in the latter case the initial transverse wave parameter K^(r\z) and the mode index n are renormalized by the 
gradient factor of (1 -f S^)~^/^. This enables us to immediately write down the dispersion equations for a DDR with 
random inhomogeneous side boundary based upon the results given in Appendix [X] Specifically, for a rough-bounded 
DDR the oscillation spectrum is governed by the equation 



1 Jn'jklR) 



Hf(k^R) 
H^i\k^R) 



i?2 



(e-1) 



2 7,2 



(fcl) ki 



(33) 
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which has to be supplemented with additional relationships between k and fcz, the latter resulting from joining the 
fields at the end interfaces. In Eq. ([55)1 . the notations are introduced (n, k, k'^^, k±) = (n, k, fc^, k±) j \J\ + S^, the 
wave numbers fcj^ and k\^ are defined in Eq. (|A4[) . 

As the additional connection between wave parameters k and k^ in the case of a perfect cylindrical resonator 
equation (|A9p is obtained for i^^-symmetric oscillations. For iJ^-antisymmetric oscillations Eq. (jAlOp holds. One can 
easily check that in going from unperturbed wave equation (|A1[) to equation ([5]) describing the DDR with rough side 
wall the equations derived through joining EM field components at the end boundaries of the dielectric disk remain 
unchanged. 

Finally, proceeding from the above calculations we are led to conclude that with any roughness of the DDR side wall, 
which is basically restricted by the smallness condition ([7]), to obtain the oscillation spectrum one can make use of the 
relationship Eq. (j33[) mainly coincident in form with the dispersion equation for an infinitely long dielectric cylinder. 
As a supplementary condition to interconnect the lengthwise and transverse wave vector components equation (|A9[) 
or (|A10[) should be applied, depending upon iJ^-symmetry of the desired solution. The fundamental difference of 
Eq. (j33p from its perfect-cylinder counterpart is the renormalization of basic wave parameters, which is governed by 
geometric properties of the asperities. A special emphasis should be placed on the fact that it is exactly the gradient 
scattering mechanism, not the amplitude one, that most crucially affects the DDR spectrum. This particular fact 
makes itself evident in that the renormalization of wave parameters in Eq. (j33p is dictated not by the height-type 
parameters, such as, e. g., Rayleigh parameter ka or the ratio a/i?, but is mainly regulated by the mean-square slope 
of the asperities against the unperturbed resonator boundary, which is specified by parameter ct/sc- 

V. EXPERIMENTAL RESULTS AND DISCUSSION 

The main goal of experimental studies in this work was to validate our theory as regards microresonator spectra and 
the effect of random surface inhomogeneities on them. The point is that essential assumption adopted in the theory 
is that it virtually does not consider electromagnetic fields radiated from the resonator external edges into the corner 
regions labeled by numbers 4 and 4' in Fig. [5] (see Appendix [X)) . At the same time, without making quite complex 
calculations one cannot make a 'priori statements about these fields being small enough to neglect them in calculating 
the resonator spectrum. Yet another goal of the experiment was to examine our theoretical findings concerning the 
physical mechanism that adequately describes the influence of random surface inhomogeneities on microresonator 
spectral properties. 

The studies of the microresonator spectrum were performed through modeling these quite small systems in the 
millimeter wave band. For this purpose we used a quasioptic dielectric disk resonator. Physically, oscillation properties 
of this macroscopic system are identical with properties of silicon microresonators used as oscillation systems in real 
optical lasers. Our resonator was made of teflon whose permittivity is not signiflcantly far from unity (e = 2.08) and 
whose dielectric loss in the millimeter range are fairly small (tan 5 = 2.3 x lO""*). In the model DDR, whispering 
gallery oscillations were excited with the EM held concentrated at the periphery of the disk, in the narrow region close 
to its side boundary. This enabled us to use the disk core to fix it in the level position without introducing additional 
dissipative losses to the experiment. The source of WG modes was positioned close to the resonator side boundary. 
Its role was played by the waveguide antenna powered by the microwave generator. The antenna was fabricated as a 
waveguide tapered along the short wall, whose butt end was positioned near the resonator side surface. The receiving 
antenna was made identical to the source antenna and was placed at the diametrically opposite point off the disk. 

Thin dielectric bracket-bars were used as the inhomogeneities of the resonator side boundary. They were attached 
to the DDR side surface. The basic requirement imposed on the inhomogeneities was for them not to cause noticeable 
additional dissipative losses. To this end, the bracket-bars were made of the same tefion as the resonator body. The 
inhomogeneities distribution on the resonator side surface was random and varied for each different realization. In 
Fig. [51 the general view of the resonator is shown along with the exciting and the receiving antennas as well as with 
the attached teflon bracket-bars. 

In order to excite TE or TM oscillations in the resonator, two different configurations of antenna- versus-resonator 
were used. For TE oscillations the antenna magnetic field was directed along the resonator axis, z. To achieve this, 
the wide plate of the waveguide was aligned parallel to this axis. In the case of TM oscillations, along the axis z the 
electric field was directed. For this purpose the wide side of the antenna was oriented transversely. 

By changing the distance between the exciting waveguide butt end and the resonator side boundary, as well as the 
angle between them, we were able to adjust the coupling between the antenna and the resonator to optimize it. As 
the optimal coupling we accepted the one whereby the additional loss caused by the antenna was much less than the 
eigen-loss in the resonator. At the same time, it was necessary to keep the level of the coupling sufficient for spectral 
lines to be traceable. The pattern of whispering gallery EM fields in the DDR in known to be very sensitive to the 
frequency variation. Therefore the resonator-to-waveguide coupling is different for each of the modes. Based upon 
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FIG. 2: Schematic view of the experimental DDR with surface inhomogeneities: 1 — teflon 
disk; 2 — superimposed teflon bracket-bars (inhomogeneities); 3, 4 — the exciting and receiving 
waveguide antennas; d — the gap between the antenna and the DDR's side surface, which is used 
to tune the coupling of antenna with the resonator. The diameter of the DDR is 102 mm, the 
thickness is 7.6 mm. The inhomogeneities dimensions are: the length, the width, and the 
thickness 7.6 mm, 3 mm, and 2 mm, respectively. 

this, we established the optimal coupling for each of the spectral lines separately while carrying out the measurements 
in a wide frequency range. 

Spectral measurements with the model resonator were made in the on-pass regime using millimeter waveband 
standing-wave ratio meter. Since the experiments were conducted over a wide range of frequencies and for large 
number of realizations of surface inhomogeneities, all spectral measurements were rendered automatic. The signal 
from the ratio meter was sent to the computer and processed by means of the especially designed program to find 
both the frequencies and the quality factors of resonance lines. The accuracy of the measurements was 0.01% and 5% 
for the resonance frequencies and the quality factors, respectively. 

To identify the spectral lines, it was necessary to determine the value of azimuth index for each of the lines in the 
absence of inhomogeneities. To this end, the miniature rotating probe made of a thin metal plate was used, which we 
inserted into the region at the resonator disk where the electric field antinode was positioned. In so doing, the source 
was tuned to the frequency of the particular spectral line. When the probe rotated about the resonator axis, the 
signal registered by the receiver varied in time at the rate the probe passed across the regions with electric field loops. 
This enabled us to determine the desired mode index through the measurements of the signal modulation frequency. 
In Fig. [21 the results of numerical calculations of the resonator spectrum are shown along with spectral measurements 
data. As seen from Fig. [3^, the measured spectra of TE and TM modes correlate well with the calculated spectra. The 
difference between spectral line frequencies found from Eq. (j33p . including intrinsic dissipation loss in the dielectric, 
and experimentally measured frequencies is no more than 1%. The quality factors obtained in the experiment (as 
shown in Fig. [8)3) are smaller than those calculated theoretically. We attribute this to the fact that the measured 
Q-factor includes not only eigen-loss in the resonator material but also is affected by the loss resulting from coupling 
with the antenna. From the diagrams in Fig. [3] one can observe the nearly equidistant character of the spectrum (the 
frequency interval = 0.7 GHz), which is typical for open resonators with WG-type oscillations 16:], as well as the 
exponential dependence on the mode index n of the calculated and the measured Q-factors. 

Our calculations suggested, and that was experimentally confirmed, that for TM oscillations the resonator quality 
factor is significantly larger than that for TE oscillations, both of them being taken with the same azimuth index. 
This suggests that in the case of TM oscillations the DDR has the property to more efficiently retain electromagnetic 
field in its volume than the resonator with TE oscillations does. We have thus established that both the theoretical 
model of the resonator we have chosen for our calculations and the characteristic equations obtained thereupon have 
found impressive experimental confirmations. 

According to the findings of our theory, the main physical mechanism for EM field scattering by random inhomo- 
geneities of the resonator side boundaries is the gradient mechanism. The scattering arisen due to fluctuations in 
the asperity slope can be allowed for by way of modification of the cylindrical functions indices along with the wave 

numbers in the characteristic equation Eq. p3p . The modification reduces to multiplying the parameters n, k, k± and 

1/2 

fc^ by factor of x = (l + S^) < 1. In order to identify the main scattering mechanism experimentally and thus 
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Azimuth index (n) Azimuth index (n) 

FIG. 3: (Color online) The collation of theoretic and experimental data for the DDR with perfect cylindrical boundaries: the 
dependence a) of spectral line frequency and b) of the logarithm of quality factor on the azimuth index n. For TMn.i mode: 1 
— calculation, 2 — experiment; for TEn,i mode: 3 — calculation, 4 — experiment. 



to check the developed theory it was necessary to estimate the value of the parameter 5 starting from the parameters 
relevant to a particular experiment. We were governed by the following considerations. Setting the correlation length 
equal to an average distance between the centers of the attached dielectric bracket-bars, i. e., Sc = 27tR/N, where N 
is the total number of the bracket-bars, with Gaussian distributed inhomogeneities we obtain 



X 



n -1/2 



(34) 



It can be easily seen that the parameter x tends to decrease with an increasing number of bracket-bars, which leads 
to a decrease in the effective wave number k and the effective mode index n. Since the dependence of the quality 
factor on the mode index is nearly exponential, Q « e", the availability of surface inhomogeneities should result in 
the decrease in the resonator quality factor, whose origin is not an additional dissipative loss. 

According to our theory, with a small number of inhomogeneities the effect of slope-controlled scattering must 
be fairly slight (the parameter S decreases with lowering N). Therefore the value of quality factor must be weakly 
dependent on N as well. Such was indeed the case in our experiments whose results are depicted in Fig. |4l In 
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20 30 40 50 60 
Number of inhomogeneities (Af) 

FIG. 4: (Color online) The Q-factor curves versus the number of inhomogeneities on the resonator side 
surface. TM44,i oscillations: 1 — theory, 2 — experiment; TE44,i oscillations; 3 — theory, 4 — experiment. 
Calculations were carried out at cr = 3.5 mm and R — 51 mm. 



the same figure, the curves Q{N) are plotted for comparison, which are calculated from characteristic equation ([55)1 
including modification factor x f^nd the dissipation loss in the resonator material. The loss was taken into consideration 
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phenomenologically, by adding the imaginary part to permittivity e, see Eq. ([5]). With the statistical nature of our 
measurements (the averaging is done over a large number of realizations of the bracket-bar set), the agreement between 
theoretical and experimental results appears to be quite satisfactory. This suggests that, since both in our theory and 
in the numerical calculations based upon it the slope scattering mechanism alone is taken into account, the qualitative 
agreement between theory and experiment unambiguously corroborates the dominant role of this particular type of 
wave-surface scattering. 

As is also seen from Fig. 21 the influence of surface inhomogeneities on the DDR spectrum differs significantly 
for TE- and TM-type oscillations. The quality factors of TE oscillations are to a far larger extent subjected to the 
resonator boundary roughness as against the factors of TM oscillations. With a given number of inhomogeneities, the 
Q-factor of TE oscillations is noticeably smaller than that of TM oscillations. This fact, which is quite important for 
practical applications of microresonators in laser oscillation systems, was given certain attention in Ref. Yet in 
the present work it has been comprehensively substantiated. 

VI. CONCLUSION 

To summarize, we have investigated spectral properties of a dielectric disk resonator with randomly rough side 
boundary both theoretically and experimentally. It is shown that the azimuth modes of the resonator oscillations can 
be rigorously separated at whatever level of surface inhomogeneities. This allowed us to obtain the asymptotically 
exact dispersion equations which are valid over a wide range of the roughness parameters. The only requirement 
imposed on the inhomogeneities and effectively utilized in the calculations was that their mean-square height be small 
as compared with the unperturbed (i.e., non-rough) disk resonator radius. 

In deriving dispersion relations, we have shown that electromagnetic wave scattering resulting from the boundary 
roughness can be described in terms of two fundamentally different physical mechanisms, specifically, the amplitude 
(the height) scattering mechanism and the gradient (the slope) mechanism. For the first of them, the ratio of the mean 
height of the asperities to the oscillation wavelength (the Rayleigh parameter) acts as the main guiding parameter 
whereas for the gradient mechanism the mean slope of the asperities relative to the unperturbed resonator boundary 
plays the same role. Our estimations have revealed that it is exactly the gradient scattering that is of primary 
importance for the formation of rough resonator spectrum. The effect of this type of scattering can be described 
through the approximate dispersion equation that differs from the analogous equation for the ideally circular DDR 
by the gradient renormalization of the basic wave parameters, i.e,, the mode wave numbers and the azimuth index. 

Our theory aimed at describing the influence of random surface inhomogeneities on microresonator spectral prop- 
erties analytically is actually of model nature. The point is that it does not incorporate the electromagnetic fields 
radiated from the resonator into external angle sectors. To neglect these fields theoretically is quite a challenge. 
Therefore, to verify the conclusions of our approximate (to the certain extent) theory we have made experimental 
measurements on the model millimeter wave band dielectric disk resonator. Surface inhomogeneities were presented 
by teflon bracket-bars of relatively small cross-section, which were attached randomly to the resonator side bound- 
ary. For the perfect-wall resonator, i. e. the resonator with no inhomogeneities attached, our measurement data 
have demonstrated excellent agreement with the developed theory concerning both the frequency spectrum and the 
quality factors of spectral lines. Close qualitative agreement is also attained with regard to the effect produced on 
the resonator spectrum by random surface inhomogeneities placed at the side boundary. First, we have confirmed the 
theoretical predictions about the leading role of the gradient scattering mechanism in describing the effect of random 
rough boundaries. Second, we ascertained that the system of model dispersion equations we have obtained herein is 
quite effective in describing the open microresonator spectra. 

Besides, our experiments have revealed that the effect produced by the surface inhomogeneities of the DDR on the 
TE and TM oscillation spectra is fundamentally different. The TM oscillations quality factor exceeds significantly the 
analogous factor for the TE oscillations, being yet less affected by surface inhomogeneities. In our work this particular 
fact, which is essential for the production of microresonator-based lasers, has been theoretically and experimentally 
validated. 
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APPENDIX A: MODEL SPECTRUM OF IDEAL CYLINDRICAL DDR 



Proceeding from the methodology considerations, we outhne the technique for deriving dispersion equations for 
cyhndrical finite-size resonator with perfectly smooth bounding surfaces. 

Consider equation Q without potentials y^''-' and V^^*^ . Upon going to Fourier representation over angle variable 
if using the complete set of eigenf unctions n) — (27r)^^/^ exp(— mt/s), where n = 0, ±1, ±2, . . . , the equation for 
n-th angular component of the wave function becomes 



1 d d 



r or or oz^ 



*„(r,z)=0. (Al) 



and Izl ^ CX3 the radiation conditions are 



The function ^„(r, z) is sought to be finite as r ^ 0, whereas at r — 
meant to be fulfilled. 

It is difficult to find the explicit solution of Eq. (|Aip in the entire domain of variables r and z. Yet basically there 
is no need to make use of such a solution. It will suffice to obtain dispersion relations which are not strictly valid but 
satisfied with good accuracy. We will seek a desired solution in the model form, as was done, e.g., in Refs. [itL [l8l|. 
imposing the pair of basic requirements. One of the requirements is to fulfil fundamental boundary conditions at 
zero and infinite distances from the resonator center, while the other is to provide correct joining of the EM field 
components at the dielectric disk interfaces. 

By representing the solution of Eq. (|Aip as a sum of iJ^-symmetric and i^x-antisymmetric summands (hereinafter 
we will specify them by indices (s) and (a), respectively) we will seek the Debye potentials in the form given below. 



Uil\r,z) 



J„(fclr)cosfc2Z , 
B';y'^Hl^\kx_r) cos k,z , 
U^^{r,z)^Cf!^ljniklr)eM~^z\z\) , 



(A2a) 
(A2b) 
(A2c) 



V!^l\r, z) = J„(fcir) sinfc.z , (A3a) 

V^^ (r, z) = Bl^^'^Hi^^ {k^r) sin fc,z , (A3b) 
T/i|)(r,z) = Cr(^)j„(fclr)exp(~>^.|z|)sgnz . (A3c) 

In Eqs. (jA2p and (jASp . the subscripts 1,2,3 correspond to the regions (see Fig. [5]) which are labeled with corresponding 
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FIG. 5: The areas on the half-plane (r, z) between which the joining of EM fields is performed. 
The shaded region is filled up with the dielectric. 



numbers. The wave parameters k±, k'^ and Xz are given by 



u2/', l\ l„2 



{k\ 



ek^ ~ kl 
k\e-l) - ki 



(A4a) 
(A4b) 
(A4c) 



By representing EM field components in terms of the potentials (jA2|) and (jA3P and by joining the in- and the out- 
field components at the resonator side boundary (i.e., between regions 1 and 2) we obtain the well-known equation 
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that couples together wave parameters k and [l6 



e Jn'{klR) 1 H^^) [k^R) 



kl UklR) k^ HL'\k^R) 



1 Jn'iklR) 1 H^^^ (k^R) 



kl UklR) k^ Hl^\k^R) 



[kifki 



(A5) 



In order to obtain yet another couphng equation for the same wave parameters we join the tangential components of 
the field at the end boundary z = H/2. This leads to the set of equations 
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Jn'{klr)k 



kz sin ^^^^A';^"'^ + exp i 
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(A6a) 



(A6b) 



(A6c) 



(A6d) 



where the number of the unknowns (to the latter we assign not only constant factors included in Eqs. (jA2p and 
(IA3P ) but also the functions J„'(fc^r) and (n/r) J„(fc^r)) exceeds the number of the equations thus obtained. The 
situation can be improved if one adds to the system (|A6p a pair of equalities resulting from joining the EM field 
normal components at the same end boundary, specifically, 



^:^^^cos^ = iQ(^)exp(-x,i//2) 
sin ^ = C,^^'') exp ( - >c,H/2) 



(A7a) 
(A7b) 



Using Eqs. (jA6p and (jA7p we arrive at a set of four coupled equations which can be naturally combined in pairs, viz. 
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(A8b) 



in two 







Since the Bessel functions do not vanish simultaneously with their derivatives, one can satisfy Eqs. 
ways. The first one is to equate the expression in parenthesis of Eq. (IA8a|) to zero, putting the constant C, 
in Eqs. (jASbp . Alternatively, the parenthesis in Eq. (IA8bp can be set equal to zero concurrently with the coefficient 
Cn^^"* in Eqs. (lASap . As can be readily seen from Eqs. (jA2p and (|A3p . the former solution corresponds to the TM-type 
oscillations whereas the latter is consistent with TE polarization. Both of these cases can be combined into one 
equation of the following form 



k 2, k-^ H 
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The lower curly brackets in Eq. (jA9[) indicate the polarization which correlates with vanishing the expression in 
the corresponding parenthesis. By carrying out the calculations similar to those given above, for i?^-antisymnictric 
solution of Eq. (jAip we arrive, instead of Eq. (|A9p . at the equality 

(^^^_fc^tan^)(^>., + :|cot^) =0. (AlO) 

TE TM 

Note that the relationships (|A5|) , (|A9p and (jAlOp were previously obtained in Ref . |l3| , although the latter two of 
them were of somewhat different form. In that paper, however, polarization of the excited EM field was not discussed 
at all. At first glance, from our derivation of Eqs. (|A9|1 and (|A10|1 it may seem that TM- and TE-type oscillations 
in the DDR of finite thickness can be perfectly separated, at least, with no regard for the boundary inhomogeneity. 
As a matter of fact, this is not the case because one must take into account the inaccuracy, i. e., the model character 
of wave solutions (|A2p and (|A3P used when deriving dispersion relations. These solutions, as well as the analogous 
ones for the antisymmetric case, are well adapted for joining the in- and the out- field components at the boundaries 
between regions 1-2 and 1-3 in Fig. [5] The fields in region 4 are not taken into consideration, which must inevitably 
result in some "overflow" of the spectrum obtained in such a way. 

It is a difficult task to correctly evaluate the possibility to neglect the fields in the 4-th region in Fig. [5] without 
resorting to rigorous calculations. Nevertheless, by now there does not exist a rigorous theory for a finite thickness 
DDR, even in the seemingly simple case where the boundaries are perfectly smooth. For this reason, in order to test 
the results we obtained by the model calculations we choose to compare the spectrum resulting from Eqs. ([33|) . (|A9|) . 
(|A10P with the one measured in the experiment. 



APPENDIX B: MODE SEPARATION IN TWO-DIMENSIONAL WAVE EQUATION WITH ARBITRARY 

SCATTERING POTENTIAL 



Consider the equation for Green function of Eq. (|13p without rendering both the physical nature of mode potentials 
and their absolute value concrete. 



Gnn'{r,r'\z) - E Unm{r\z^)Gjnn'{r,r'\z^) = -S{r - r')Snn' ■ (Bl) 



ra^n 



Along with the exact Green function, which has a matrix structure in variables n and r, we introduce for each of the 
azimuth modes the trial mode propagator, which is assumed to obey the closed equation 



G'y\r,r'\z^)^-5{r-r') (B2) 



resulting from Eq. (jBip providing that the inter-mode scattering is disregarded. We will use the notation Gi^"* ^ for 



the operator in square brackets of Eq. (IB2p . The initial equation (|Bip can then be recast as 

G'^p-^Gnn'{ry\z^) - U{t - r')<5„„' + l^nUr\z^)Gmn'{ry\z^) , (B3) 

7n^n 

or, in the equivalent integral form, as 

oo 

Gnu'{r,r'\z^) = G'^;:^\r,r'\z^)5nn' + J2 J ndriG^/\r,n\z^)Unrniri\z^)Gmn'iri,r'\z^) . (B4) 

By setting in Eq. (jB4p index n ^ n' and then re-labelling all mode indices we can write this equation as the equation for 
solely inter-mode components of the Green matrix ||Gmn||, the intra-mode ones being thought of as known functions. 



Gmn{r,r'\z^) ~ E / ridriGl^\r,ri\z^)Ujnk{ri\z^)Gkn{ri,r'\z^) = / ridriGl^\r,ri\z^)Umn{ri\z)Gnn{ri,r'\z^) . 

(B5) 



17 



At this stage we introduce three operators, U and R, which are assumed to act in the reduced coordinate- 

mode space M„ consisting of the half-axis r ^ and the entire set of mode indices except for particular index n. The 
operators are specified by their matrix elements 

(r,n|^(^)|r',m) = G'^^ {t,t'\z^)6,,^ , (B6a) 

{r,n\U\r',m) ^ U„m{r\Zf;)-6{r - r') , (B6b) 

r 

{r,n\R\r',m) = Gl^\r,r'\z)Lin„r{r'\z^) . (B6c) 

Equation (|B5P can now be recast as the matrix equality 

{r,m\{l~R)g\r',n) ^ {r,m\RPng\r',n) , (B7) 

or, equivalently, in the general operator form as Pm(l ~ R}GPn ~ Pm^PnGPn- Here, P„ is the projection operator 
whose action reduces to assigning the value n to the nearest mode index of an arbitrary operator standing adjacent to 
it, either to the left or right. Multiplying both sides of the operator equality thus obtained by the operator (l — R) 
whose regularity was substantiated in Ref. [20] we arrive at the integral relation between the non-diagonal and diagonal 
mode matrix elements of the Green function, viz. 

Grnn{r,r'\zg) = {r,m\{l~Ry^RPj\r',n) = {KrnnG„n){r,r'\z^) . (B8) 

Setting in Eq. (|Bip index n' = n and substituting the inter-mode propagators in the form (|B8|) we eventually obtain 
the closed equation for intra- mode Green function G„„(r, r'|zg), 

Gnn{r,r'\z^)^-6(r-r') . (B9) 
r 

Here, 7^ is the operator that accounts for the inter-mode scattering (it is just the T-matrix well-known in quantum 
scattering theory [IEHIDj which has the form 

%^ PnU{t-R)''RPn . (BlO) 

Trial Green function Gli^\r,r'\z^) entering operator potential (jBlOp through matrix elements (jB6cP obeys equation 
(IB2p and the same boundary conditions just like the desired mode propagator G„„(r, r'\z^). 

Equation (|B9|) along with Eq. (|B8|) determines the entire Green function of wave equation (fT3|) and serves as a basis 
for deriving the system of uncoupled equations for azimuth components of the sought-for wave function. 



1 d d / I N " 

-^^^ + Klir\z^) - -. 
r or or r 
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